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We have recently shown [Blunt, MO, et al. (2008) Random tiling 
and topological defects in a two-dimensional molecular network. Sci- 
ence, 322:1077-1081] that p-terphenyl-3,5,3',5'-tetracarboxylic acid 
adsorbed on graphite self-assembles into a two-dimensional rhom- 
bus random tiling. This tiling is close to ideal, displaying long range 
correlations punctuated by sparse localised tiling defects. In this 
paper we explore the analogy between dynamic arrest in this type 
of random tilings and that of structural glasses. We show that the 
structural relaxation of these systems is via the propagation-reaction 
of tiling defects, giving rise to dynamic heterogeneity. We study the 
scaling properties of the dynamics, and discuss connections with ki- 
netically constrained models of glasses. 

random tiling | glass transition | dynamic heterogeneity 

Abbreviations: KCM, kinetically constrained model; STM, scanning tunneling micro- 
scope; TPTC, p-terphenyl-3,5,3',5'-tetracarboxylic acid 



Introduction 

Figure 1 shows a molecular network formed by organic 
molecules adsorbed from solution onto a graphite substrate 
1 . The molecule, p-terphenyl-3,5,3',5'-tetracarboxylic acid, 
or TPTC (see Fig. 1A), binds to other TPTC molecules on the 
substrate adopting one of three possible orientations. Each 
molecule can then be mapped to a rhombus tile, see Fig. IB, 
where the colours red, green and blue indicate the molecular 
orientation. Neighbouring molecules (tiles) can bind to neigh- 
bours in a parallel or "arrowhead" configuration, equivalent to 
junctions between tiles of the same or different colour, respec- 
tively. Fig. 1C shows a scanning tunneling microscope (STM) 
image of the resulting molecular network of adsorbed TPTC, 
and Fig. ID the corresponding rhombus tiling [I] where each 
molecule is represented by a tile. 

The molecular networks studied in [1] are close to "per- 
fect" rhombus tilings (or dimer coverings of the honeycomb 
lattice) [2j [3, 4 , in the sense that they contain rather few 
tiling defects, typically less than one defect per 300 adsorbed 
molecules pp. In Figs. 1C and ID one such tiling defect is iden- 
tified. They are also entropically stabilized "random tilings" 
displaying algebraic spatial correlations [T], characteristic of 
a critical, or Coulomb, phase [2[3l|4]. The structures such 
as those of Fig. 1 are close to dynamically arrested at room 
temperature pQ. The interaction energy between two neigh- 
bouring molecules is several times ksT [I], so once a tiling is 
formed tile removal is highly suppressed, and structural relax- 
ation is slow. Tile rearrangements mediated by propagation 
of defects have been observed experimentally, but so far on 
timescales of seconds for each event ]T]. This combination of 
an amorphous structure, albeit with critical spatial correla- 
tions, and very slow relaxation suggests an analogy between 
this kind of random tilings and structural glass formers [6]. 

The aim of this paper is to discuss this analogy. For sim- 
plicity we focus on the case where all tile-tile interactions are 
equal, since the dynamics for small bias [7] is qualitatively 
the same, and on the long-time dynamical regime; the initial 
growth dynamics is discussed in a separate paper [TO]. We 
show, by means of numerical simulations, that the low tem- 
perature dynamics of a rhombus tiling where the number of 
tiles is not conserved displays some of the features observed 



in the dynamics of structural glass formers, in particular dy- 
namic heterogeneity. Relaxation in these random tilings is 
facilitated by tiling defects, a mechanism similar to that of 
kinetically constrained models of glasses. We will also discuss 
this connection. 
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Fig. 1. (A) Molecular ball and stick diagram of p-terphenyl-3,5,3,5-tetracarboxylic 
acid (TPTC). (B) Example of an arrangement of TPTC molecules linked via hydro- 
gen bonds when adsorbed on substrate, and rhombus tile representation; the tiles 
are coloured according to the three possible orientations of the molecule. (C and 
D) Mapping to a Rhombus Tiling: (C) shows an STM image of a typical area of 
TPTC network adsorbed on graphite; the backbones of the TPTC molecules appear 
as bright rods in the image. The corresponding rhombus tiling is shown in (D). The 
molecular network is a rhombus tiling of the plane or, equivalently, a dimer covering 
of the honeycomb lattice. The white circle shows the position of a tiling defect. (E 
to L) Example of motion of tiling defects: The leftmost defect (upward pointing 
triangle) effectively makes two hopping steps, between (E) and (G), and (G) and (I). 
This motion is mediated by the desorption (F and H) and re-adsorption of a tile (G 
and I). Between (I) and (K) the rightmost defect (downward pointing triangle) makes 
a step to the left, which brings it into contact with the leftmost defect. The two 
annihilate with the adsorption of the last tile (L). (See Ref. 1 for details.) 
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Fig. 2. (A) Relaxation of the concentration of defects c(t) as a function of time, 
starting from an empty lattice, c(0) = 1, for various temperatures T . (B) Same as 
(A), but time rescaled by the defect effective hopping rate V ~ e — The curves 
collapse at long times in this representation. The dotted line indicates the power-law 
decay (Tt)~ ^'^ expected from diffusion-pair-annihilation, A-\-B — > 0, in d = 2. 
The observed behaviour is closer to c(t) ^ (Tt)~ as indicated by the dashed 
line. 



Model 

We simulate a dimer covering of the honeycomb lattice, which 
is equivalent to a rhombus tiling of the plane [2] [3]. That is, 
each rhombus tile is composed of an upward and a downward 
pointing triangle face-to-face; these triangles are centred at 
the sites of a honeycomb lattice, each in a different sublat- 
tice. The dynamics we consider is one where the only possible 
moves are the adsorption of a tile on the lattice, if the two 
sites it would occupy are empty, or the desorption of a tile. 
The number or tiles (and of tiling defects) is therefore not 
conserved. This resembles the experimental situation where 
molecules are exchanged between substrate and solution pQ. 
We only consider single-tile moves [12] . We set the binding 
energy J between neighbouring tiles to J = 1. For low tem- 
peratures, T < 1 (in units of J, and where fca = 1), we sim- 
ulate the dynamics using a version of Borz-Kalos-Lebowitz, 
or continuous-time, Monte Carlo [141 I15j which is particu- 
larly efficient for this problem. We simulate systems of sizes 
N = 200 x 200 to N = 10 3 x 10 3 at all temperatures. 

At low temperatures, when the density of tiles is high, 
desorption of tiles is rare. The energy barrier to remove a tile 
surrounded by four neighbouring tiles is AE = 4, and the rate 
for that transition is suppressed by a factor of e~ i ^ T . A more 
likely transition is the removal of a tile neighbouring a defect, 
as the rate for this process scales as e _3//T . This will give rise 
to the effective propagation of tiling defects, as sketched in 
Figs. 1E-1L. In Fig. IE there are two tiling defects of opposite 
"charge" [4]. Fig. IF shows the desorption of a tile next to 
the leftmost defect (a process of rate oc e~ 3 ^ T ) and Fig. 1G 
the subsequent adsorption of another tile (a process of rate 
O(l), as it is energetically favourable). The net effect is the 
hopping of the defect by one (sub)lattice site. The sequence 
1G-1I shows a second such step. The effective hopping rate is 
therefore F oc e _3//T . 

Defects of opposite charge can effectively react with each 
other. This is sketched in Figs. 1I-IL. In this case the right- 
most defect hops one step to the left. On meeting the opposite 
defect a gap large enough for a tile is formed, and the two de- 
fects are "annihilated" by the adsorption of a tile, at rate 
O(l) as this is energetically favourable. Of course, this pro- 
cess is reversible, and two opposing defects can be "created" 
by desorption of a molecular tile, at a rate oc e~ 4 / T . The ef- 
fective dynamics of defects, therefore, resembles a reversible 
A + B <-> reaction-diffusion process |16|, although, as we will 



see below, it is not clear that defect propagation is actually 
diffusive or that defect interaction can be approximated as 
being local (see also [T71 \TE\ ). 



Results 

Figure 2A shows the evolution in time of the concentration 
of defects, c(t), starting from an empty lattice at time zero, 
c(0) = 1, at various temperatures T. After a short initial 
transient of fast, temperature independent, tile adsorption, 
the system enters a regime of activated dynamics: most de- 
fects are isolated, and energy barriers need to be crossed for 
the dynamics to progress. The dynamics becomes increas- 
ingly sluggish with decreasing temperature, and once times 
are long enough for defect motion to take place relaxation en- 
ters a scaling regime. Fig. 2B shows that the rate limiting 
step is defect hopping: the long time data collapses if time 
is rescaled by the defect hopping rate, t — > Ft. The defect 
concentration decays as c(t) ~ (Tt)~ a , with a « 3/4. This 
exponent is somewhat different from the exponent a = 1/2 of 
two-species diffusion-annihilation, A + B — > 0, in dimension 
two [161 119| . This could be an indication that defect prop- 
agation is non-diffusive, although an exponent of a = 3/4 
can also be explained by initial state fluctuations in the tiling 
case which differ from those of the standard A + B — * prob- 



lem [20]. Eventually, for times t S> 



„4/T 



the reverse process 



— > A + B becomes accessible, and the concentration relaxes 
to its equilibrium value c(t) — + c cq . 
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Fig. 3. (A) Equilibrium concentration of defects, c C q, as a function of tem- 
perature T. The straight line corresponds to the fit c C q = e — "V^. (B) 
Equilibrium height correlations at various temperatures. The main panel shows 
{[h(r) — h(0)] 2 ) as a function of distance r. As the defect concentration de- 
creases with decreasing temperature, the curves approach the ideal tiling behaviour 
{[h(r) — h(0)} 2 ) = 9/-7T 2 In (r). The Inset shows the correlation function 
( e ikAh(r)y where A/i(r) = h(r) - h(0), for k = tt/5; the correlation 
behaves in a similar way for other choices of h. The dashed line is the power law 
behaviour for the ideal tiling, r~ 9k / 27r = r~ 9 / 50 for this choice of k. Once 
again, the lower the temperature, the longer the algebraic regime. For one of the 
temperatures, T = 0.4, we show that < e * fcA M r )} = e ~k 2 /2([h(r)-h(0)f) _ 
as expected from a Gaussian form of the free-energy for the height. (C) Scheme for 
obtaining the height representation of a tiling. A displacement along a tile edge leads 
to an increase in height by +1 or —1 as shown (cf. Fig. 1). 



2 | www.pnas.org 



Footline Author 



The equilibrium properties of the tilings are shown in Fig. 
3. The temperature dependence of the equilibrium defect 
concentration is given by c cq ~ e _3/ ^ T , see Fig. 3A. This is 
what one would obtain for a non-interacting gas of defects 
on the lattice with an energy cost of E = 3 per defect. Fig. 
3B shows the spatial correlations. A rhombus tiling can be 
mapped to a height field on the triangular lattice [3]: the 
height h changes by ±1 unit when traversing the edges be- 
tween tiles, according to the prescription of Fig. 3C. The 
main panel of Fig. 3B shows the height-height correlation 
function, (\h(r) - h(0)] 2 }, as a function of distance r (along 
lattice directions), for various temperatures [23) . At low tem- 
peratures this correlation approaches the ideal tiling limit 
{[h(r) — h(0)] 2 ) = 9/7r 2 ln(r), corresponding to a Gaussian 
free-energy F — J d 2 x(K/2)\ X7h(x)\ 2 for a continuous height 
field U], with elastic constant K = n/9 [3]. For finite T the 
logarithmic behaviour is over a finite distance due to the pres- 
ence of tiling defects [24] . An alternative correlation function, 
( e ife[fe(r)-Ma)]) [22] , is shown in the Inset to Fig. 3B, for the 
specific value k = 7r/5 of the "height space" reciprocal vector 
|22| (the behaviour is similar for other choices of k). At low 
temperatures the function becomes algebraic indicating long 
range correlations. The decay exponent is close to 9k 2 /2ir 2 , 
as expected from the Gaussian form of the free-energy |22| . 
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Fig. 4. Dynamic heterogeneity in random rhombus tilings. The panels show the 
persistence field Pi(t) of the local autocorrelation function Ciit) (see text) at var- 
ious times t, for representative equilibrium trajectories, at two different temperatures 
T. Black indicates Pi(t) = 1, and white Pi(t) = 0. The average relaxation 
time is T, see Fig. 5. Relaxation is clearly heterogeneous. The size of dynamic 
heterogeneity grows with decreasing temperature. 



As described above, at low temperatures structural re- 
arrangements are most likely in the neighbourhood of tiling 
defects. This gives rise to heterogeneous relaxation, as il- 
lustrated in Fig. 4. Here we plot the local autocorrelation 



function d(t) 



where m stands for the state 



of site i in the lattice, say m = 0,1,2,3 for empty, or 
occupied by a red, green or blue tile, respectively. More 
precisely, Fig. 4 shows the corresponding persistence field, 
~ rit'=oC»(*) : ^ s ite i has never relaxed up to time 
t then Pi(t) = 1, and Pi(t) = otherwise. The different pan- 
els show how relaxation is distributed in space at different 
times. Clearly, the system relaxes in a heterogeneous, spa- 



tially correlated manner. Fig. 4 also suggests that the size 
of these spatial dynamic correlations grows with decreasing 
temperature. This is very similar to dynamic heterogeneity 
in structural glass formers |25| . 

Figure 5 quantifies equilibrium relaxation and dy- 
namic heterogeneity. Fig. 5A shows the average (con- 
nected and normalised) autocorrelation function C(t) 



(Ci(oo)) 



+ (1 



where A 

the autocorrelation function decays 



((C i (t))-A)/(1-A), 
c e q) 2 /3. As expected, 
more slowly the lower the temperature. The characteristic 
timescale for relaxation, r, obtained from these correlations is 
approximately r = roe A//T , with A w 6.6, see Inset to Fig. 5B. 
Relaxation times thus increase with decreasing temperature 
following an Arrhenius law. In the context of the glass tran- 
sition this is often termed "strong" glass forming behaviour 
[6]. Moreover, the autocorrelations are close to exponential, 
rather than stretched exponential [B]. This could mean that 
relaxation is not collective, but Fig. 4 suggests otherwise. The 
exponent A appears non-trivial: we have that A < 7, the 
energy barrier to remove two neighbouring and parallel tiles 
(which would allow the adsorption of a distinct tile in the 
space created); this is the lowest energy barrier to purely local 
relaxation of the autocorrelation. This indicates that relax- 
ation is achieved more effectively through defect propagation, 
a collective mechanism. Interestingly, we also have A > 6, 
which is the value one would expect for diffusing defects, al- 
though this may simply be due to logarithmic corrections to 
diffusion in dimension two 1261. 
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Fig. 5. (A) Average equilibrium autocorrelation function C(t) at various temper- 
atures T. (B) Four-point susceptibility \4 {t) at the same temperatures as (A). The 

four-point susceptibility has a maximum at t ~ T. The peak value X4 P ° a ^^ ~ t\ 
with A ftrf 0.4, see dotted line. The growth towards the peak follows the power law 
X4 ^ , with fi Pi 1.6. The Inset to (B) shows that the relaxation timescale for 
the autocorrelation function follows an Arrhenius law, r oc e^'^, with A 6.6. 



In Fig. 5B we show the "four-point" susceptibility, X4 W = 
N(l - A)' 2 [( J2a Ci(t)Cj(t)) - C 2 (t)] , which measures 
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sample to sample fluctuations in the correlator C(t). This 
is an observable often used to quantify dynamic heterogeneity 
|25l 1271 [28] . X4 is non-monotonic in time, peaking at times 
close to the structural relaxation time, t ss r, where dynamic 
heterogeneity is most prominent. The size of the peak of \4 in- 
creases with decreasing temperature, indicating that dynam- 
ical fluctuations are larger at lower T. The four-point sus- 
ceptibility displays dynamic scaling, but again this is slightly 
different from what one expects from diffusing excitations in 
two dimensions [27]: the peak value, X4° ak = Xi( T )j scales 
as a power of the relaxation time, x^T^ ~ rA > with A « 0.4 
(rather than X = 1); and the growth of X4 towards the peak 
goes as Xi{i) ~ i M , with n m 1.6 (rather than n = 2). 

Discussion 

We have shown that random tilings, of the kind correspond- 
ing to the experimental system of 1 , display features com- 
monly associated with glass forming systems, most notably a 
pronounced slowdown at low temperatures and accompanying 
heterogenous relaxational dynamics. In these nearly dynami- 
cally arrested tilings, structural relaxation occurs through the 
propagation of rare localised tiling defects. This is an exam- 
ple of the mechanism of dynamic facilitation [29[ |3"D] , whereby 
defects "facilitate" molecular rearrangements in their imme- 
diate vicinity. The fact that defects are scarce at low T, and 
that their motion is activated, leads to the observed slow- 
down, and to fluctuation-dominated, heterogeneous dynam- 
ics [30j . Relaxation by means of defect propagation-reaction 
is the hallmark of kinetically constrained models (KCMs) of 
glasses [11] . The dynamics of the random tiling systems stud- 



ied here is very close to that of these idealised models. For the 
simple rhombus tiling, when the number of tiles is not con- 
served, we have found Arrhenius timescales and exponential 
relaxation. This is similar to the simpler KCMs, such as the 
Fredrickson- Andersen model [11] in dimensions two or more 
[31] . The observed dynamic scaling properties may suggest 
however that defect dynamics is not simply that of diffusion- 
annihilation-creation, but this requires further exploration. 

It is possible that more complex random tiling sys- 
tems, such as those giving rise to quasicrystals [4 , display 
even richer slow dynamics, in particular, super- Arrhenius 
timescales at low temperatures and stretched relaxation func- 
tions [32| . Random tilings also offer a further testing ground 
for theories of the glass transition. Specifically, it would be of 
interest to see whether their slow dynamics can be explained 
using the "mosaic" perspective on glasses [361 137]. Here we 
have a system with a finite configurational entropy density 
given by all the possible local tiling arrangements, the central 
tenet of the mosaic approach [36] [37]. The mismatch between 
mosaics however does not produce extended interfaces, but 
rather localised defects, which as shown above make the dy- 
namics close to that of KCMs. Can one therefore explain the 
structural relaxation of these systems via the entropic droplet 
picture of Refs. [36 | I37]? And could this highlight connections 
between the mosaic approach and that based on dynamic fa- 
cilitation and KCMs [38]? We hope to address these and other 
questions in future work. 
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